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Abstract 

A measurement of the double-differential ir ± production cross-section in proton-carbon, proton- 
copper and proton-tin collisions in the range of pion momentum 100 MeV/c < p < 800 MeV/c and 
angle 0.35 rad < 9 < 2.15 rad is presented. The data were taken with the HARP detector in the 
T9 beam line of the CERN PS. The pions were produced by proton beams in a momentum range 
from 3 GeV/c to 12 GeV/c hitting a target with a thickness of 5% of a nuclear interaction length. 
The tracking and identification of the produced particles was done using a small-radius cylindrical 
time projection chamber (TPC) placed in a solenoidal magnet. An elaborate system of detectors in 
the beam line ensured the identification of the incident particles. Results are shown for the double- 
differential cross-sections d 2 a/dpd8 at four incident proton beam momenta (3 GeV/c, 5 GeV/c, 
8 GeV/c and 12 GeV/c). 
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1 Introduction 



The HARP experiment [T] makes use of a large-acceptance spectrometer for a systematic study of hadron 
production on a large range of target nuclei for beam momenta from 1.5 GeV/c to 15 GeV/c. The 
main motivations are the measurement of pion yields for a quantitative design of the proton driver of 
a future neutrino factory [2J, a substantial improvement of the calculation of the atmospheric neutrino 
flux O |H O [6] [7] and the measurement of particle yields as input for the flux calculation of accelerator 
neutrino experiments, such as K2K [8] [9], MiniBooNE [10] and SciBooNE [FT] . 

The measurement of the double-differential cross-section, d 2 a n /&p&6 for tt^ production by protons of 
3 GeV/c, 5 GeV/c, 8 GeV/c and 12 GeV/c momentum impinging on a thin carbon, copper or tin target 
of 5% nuclear interaction length are presented. 

Especially for carbon targets it is interesting to measure pion production cross-sections in the framework 
of the HARP measurement programme for neutrino flux calculations. Carbon targets are frequently used 
as hadron production targets in neutrino beam lines. In addition, measurements on carbon can be used 
to predict pion production off nitrogen and oxygen nuclei without a large extrapolation in the production 
models. The knowledge of the latter production cross-sections are needed to model atmospheric muon and 
neutrino fluxes. Owing to the relatively low incoming beam momenta the data are especially interesting 
for the calculation of hadron production in secondary interactions in extended production targets and in 
atmospheric flux calculations. The comparison of the measurements on copper and tin targets with the 
carbon data in this paper and with the tantalum data obtained with the same apparatus described in 
Ref. [12] can be used to check the dependence on the atomic number A in hadron production models. 
Copper and tin are interesting target materials as their atomic numbers are midway between light target 
materials, such as Be, Al and C (used in targets for conventional neutrino beams) and heavy targets such 
as Ta (relevant for the optimization of neutrino factory designes). 

Data were taken in the T9 beam of the CERN PS. For this analysis, about 1,159,000 (1,066,000 and 
1,284,000) incoming protons were selected which gave an interaction trigger in the Large Angle spec- 
trometer, resulting in 235,000 (209,500 and 243,400) well-reconstructed secondary pion tracks for the 
carbon (copper and tin) target. The different settings have been taken within a short running period so 
that in their comparison detector variations arc minimized. 

The analysis proceeds by selecting tracks in the Time Projection Chamber (TPC) in events with incident 
beam protons. Momentum and polar angle measurements and particle identification are based on the 
measurements of track position and energy deposition in the TPC. An unfolding method is used to 
correct for experimental resolution, efficiency and acceptance and to obtain the double-differential pion 
production cross-sections, with a full error evaluation. A comparison with available data is presented. 
The analysis follows the same methods as the ones used for the determination of ^ production cross- 
sections by protons on a tantalum target described in Ref. [T2]. We refer to Ref. 42] for a detailed 
description of the analysis, only the main points and differences with respect to the latter are described 
here. 



2 Experimental apparatus 

The HARP detector is shown in Fig. [1] The forward spectrometer is built around a dipole magnet for 
momentum analysis, with large planar drift chambers (NDC) [13J for particle tracking and a time-of-flight 
wall (TOFW) [14], a threshold Cherenkov detector (CHE) and an electromagnetic calorimeter (ECAL) 
used for particle identification. The forward spectrometer covers an acceptance for tracks originating 
from the target with polar angles up to 250 mrad. This is well matched to the angular range of interest 
for the measurement of hadron production to calculate the properties of conventional accelerator neutrino 
beams [TBI HZ]- In the large-angle region a cylindrical TPC with a radius of 408 mm is positioned in a 
solenoidal magnet with a field of 0.7 T. The target is inserted into the inner field cage of the TPC. The 
TPC is used for tracking, momentum determination and the measurement of the energy deposition AE /dx 
for particle identification [18]. A set of resistive plate chambers (RPC) form a barrel inside the solenoid 
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Drift Chambers 




Figure 1: Schematic layout of the HARP detector. The convention for the coordinate system is shown in 
the lower-right corner. The three most downstream (unlabelled) drift chamber modules are only partly 
equipped with electronics and are not used for tracking. 

around the TPC to measure the arrival time of the secondary particles [19J . Beam instrumentation 
provides identification of the incoming particle, the determination of the time when it hits the target, 
and the impact point and direction of the beam particle on the target. Several trigger detectors are 
installed to select events with an interaction and to define the normalization. 

In addition to the data taken with the thin carbon, copper and tin targets of 5% nuclear interaction length 
(Ai), runs were also taken with an empty target holder, a thin 2% Ai target and a thick 100% Aj target. 
Data taken with a liquid hydrogen target at 3 GeV/c, 5 GeV/c and 8 GeV/c incident beam momentum 
together with cosmic-ray data were used to provide an absolute calibration of the efficiency, momentum 
scale and resolution of the detector. Moreover, tracks produced in runs with Pb and Ta targets in the 
same period and with the same beam settings were used for the calibration of the detector, verification 
of the event reconstruction and analysis procedures (see Ref. p~2] for further details). 

The momentum of the T9 beam is known with a precision of the order of 1% |20j . The absolute normaliza- 
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tion of the number of incident protons was performed using 300 000, 240 000 and 280 000 'incident-proton' 
triggers for the carbon, copper and tin data, respectively. These are triggers where the same selection on 
the beam particle was applied but no selection on the interaction was performed. The rate of this trigger 
was down-scaled by a factor 64. A cross-check of the absolute normalization was provided by counting 
tracks in the forward spectrometer. 

A detailed description of the HARP apparatus is given in Ref. [15]. In this analysis the detector com- 
ponents of the large-angle spectrometer and the beam instrumentation are employed and are briefly 
summarized in the following. 

A set of four multi-wire proportional chambers (MWPCs) measures the position and direction of the 
incoming beam particles with an accuracy of «1 mm in position and w0.2 mrad in angle per projection. 
At low momenta the precision of the prediction at the target is limited by multiple scattering. A beam 
time-of-fiight system (BTOF) measures time difference of particles over a 21.4 m path-length. It is made of 
two identical scintillation hodoscopes, TOFA and TOFB (originally built for the NA52 experiment [2Tjl. 
which, together with a small target-defining scintillator (TDS), also used for the trigger, provide particle 
identification at low energies. This provides separation of pions, kaons and protons up to 5 GeV/c and 
determines the initial time at the interaction vertex (to). The timing resolution of the combined BTOF 
system is about 70 ps. A system of two N2-filled Cherenkov detectors (BCA and BCB) is used to tag 
electrons at low energies and pions at higher energies. The electron and pion tagging efficiency is found 
to be close to 100%. The proton fraction in the incoming beam varies from 35% at 3 GeV/c to 92% at 
12 GeV/c. The length of the accelerator spill is 400 ms with a typical intensity of 15 000 beam particles 
per spill. The average number of events recorded by the data acquisition ranges from 300 to 350 per spill 
for the four different beam momenta. 

The target is placed inside the inner field cage (IFC) of the TPC such that, in addition to particles 
produced in the forward direction, backward-going tracks can be measured. All three targets have a 
nominal thickness of 5% Ai and a cylindrical shape with a nominal diameter of 30 mm. The 99.95% pure 
carbon target used for the measurement described here has a thickness of 18.94 mm with a variation of 
±0.02 mm. Its density was measured to be 1.88 g/cm 3 . The copper target has a purity of 99.99% with 
a thickness of 7.52 mm with a variation of ±0.01 mm and a density of 8.92 g/cm 3 . The tin target has a 
purity of 99.99% with a thickness of 11.04 mm with a variation of ±0.04 mm and a density of 7.29 g/cm 3 . 
A set of trigger detectors completes the beam instrumentation: a thin scintillator slab covering the 
full aperture of the last quadrupole magnet in the beam line to start the trigger logic decision (BS); a 
small scintillator disk, TDS, positioned upstream of the target to ensure that only particles hitting the 
target cause a trigger; and 'halo' counters (scintillators with a hole to let the beam particles pass) to veto 
particles too far away from the beam axis. A cylindrical detector (inner trigger cylinder, ITC) made of six 
layers of 1 mm thick scintillating fibres is positioned inside the inner field cage of the TPC and surrounds 
the target. It provides full coverage of the acceptance of the TPC. The large-angle spectrometer consists 
of a TPC and a set of RPC detectors inside the solenoidal magnet. The TPC detector was designed to 
measure and identify tracks in the angular region from 0.25 rad to 2.5 rad from the beam axis. Charged 
particle identification (PID) can be achieved by measuring the ionization per unit length in the gas 
(dE/dx) as a function of the total momentum of the particle. Additional PID can be performed through 
a time-of-ffight measurement with the RPCs. 

In the present analysis, the TPC provides the measurement for the pattern recognition to find the particle 
tracks, and to measure their momentum through the curvature of their trajectory. It also provides PID 
using the measurement of energy deposition. The RPC system is used in this analysis to provide a 
calibration of the PID capabilities of the TPC. 

In addition to the usual need for calibration of the detector, a number of hardware shortfalls, discovered 
mainly after the end of data-taking, had to be overcome to use the TPC data reliably in the analysis. The 
TPC is affected by a relatively large number of dead or noisy pads and static and dynamic distortions of 
the reconstructed trajectories. Static distortions are caused by the inhomogeneity of the electric field, due 
to an accidental mismatch between the inner and outer field cage (powered by two distinct HV supplies) 
and other sources. Dynamic distortions are caused instead by the build-up of ion-charge density in 
the drift volume during the 400 ms long beam spill. All these effects were fully studied and available 
corrections are described in detail in Ref. [12] . While methods to correct the dynamic distortions of the 
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TPC tracks are being implemented, a pragmatic approach has been followed in the present analysis. Only 
the events corresponding to the early part of the spill, where the effects of the dynamic distortions are 
still small, are used 0. The time interval between spills is large enough to drain all charges in the TPC 
related to the effect of the beam. The combined effect of the distortions on the kinematic quantities used 
in the analysis has been studied in detail and only that part of the data for which the systematic errors 
can be assessed with physical benchmarks was used, as explained in [TJ]. More than 40% of the recorded 
data can be used on average in the current analysis. 

The absolute scale of the momentum determination is determined using elastic scattering data off a 
hydrogen target. The angle of the forward scattered particle (pion or proton) is used to give an absolute 
prediction for the momentum of the recoil proton. This prediction is compared with the measurement 
in the TPC. To study the stability of this measurement protons are selected in a narrow band with a 
relatively large dE/dx where the dE/dx depends strongly on momentum. The average momentum for 
the protons selected in this band remains stable within 3% as a function of time-in-spill over the part of 
the spill used for the analysis. 



3 Data selection and analysis 

The beam of positive particles used for this measurement contains mainly positrons, pions and protons, 
with small components of kaons, deuterons and heavier ions. Its composition depends on the selected 
beam momentum. The analysis proceeds by first selecting a beam proton hitting the target, not accom- 
panied by other tracks. Then an event is required to be triggered by the ITC in order to be retained. 
After the event selection the sample of tracks to be used for analysis is defined. Tracks are only considered 
if they contain at least twelve space points out of a maximum of twenty. This cut is applied to ensure 
a good measurement of the track parameters and of the dE/dx. Furthermore, a quality requirement is 
applied on the fit to the helix. The latter requirement introduces a very small loss of efficiency. For 
tracks satisfying these conditions, a cut is made on d' , the distance of closest approach to the extrap- 
olated trajectory of the incoming beam particle in the plane perpendicular to the beam direction and 
z' , the z-coordinate where the distance of the secondary track and the beam track is minimal. Finally, 
only tracks with momentum in the range between 100 MeV/c and 800 MeV/c are accepted. In addition, 
particles with transverse momentum below 55 MeV/c are removed. 

Table [1] shows the number of events and tracks at various stages of the selection. The total number 
of events taken by the data acquisition ( "Total DAQ events" ) includes triggers of all types as well as 
calibration events; the number of "Protons on target" represents the count of the incoming beam trigger 
after off-line selection of accepted protons multiplied by the down-scale factor 64. The number of accepted 
events for this analysis ("Accepted protons with LAI (Large Angle Interaction)") is obtained using the 
same selection of incoming protons in coincidence with a trigger in the ITC. The large difference between 
the rows "Total DAQ events" and "Accepted protons with LAI" is due to the relatively large fraction 
of pions in the beam and to the larger number of triggers taken for the measurements with the forward 
spectrometer. These data will be the subject of other publications. The line "Maximum iV GV t" refers to the 
last number of events N cvt in spill used to avoid dynamic distortion corrections, with the corresponding 
number of interaction triggers used in the analysis ( "LAI in accepted spill part" ) and the fraction of 
the data used given under "Fraction of triggers used" . The lines "Accepted momentum determination" 
and "In kinematic region and originating from target" give the number of tracks passing the momentum 
fit quality requirements and the selection of tracks originating in the target region. Finally, the rows 
"Negative particles" , "Positive particles " , u ir~ selected with PID" and "tt + selected with PID" show the 
number of accepted tracks with negative and positive charge and the ones passing in addition the pion 
PID criteria, respectively. 

To give an impression of the complexity of the events, one can define an 'average multiplicity' as the ratio 
of the number of tracks with at least twelve hits in the TPC (regardless of their momentum, angle or 
spatial position) and the number of events accepted by the selection criteria with at least one such track. 
With this definition, the average multiplicity is 2.2, 2.6, 3.1 and 3.4 for the 3 GeV/c, 5 GeV/c, 8 GeV/c 

1 this translates into a cut on the maximum number of events (N ev t) to be retained 



7 



Table 1: Total number of events and tracks used in the carbon, copper and tin 5% Ai target data sets, 
and the number of protons on target as calculated from the pre-scaled trigger count. 



Data set 




3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


Total DAQ events 


(C) 


1304255 


2648351 


1878590 


1875610 




(Cu) 


992549 


2166883 


2599056 


748123 




(Sn) 


1636933 


2827930 


2780036 


950582 


Protons on target 


(C) 


1107456 


4872896 


6143552 


7393024 


(selected min. bias x 64) 


(Cu) 


971840 


3626048 


7606272 


2990656 




(Sn) 


1379008 


4598848 


8260544 


3842112 


Acc. protons with LAI 


(C) 


56712 


255922 


337150 


509713 




(Cu) 


59873 


237894 


541852 


226250 




(Sn) 


83549 


304949 


600581 


295053 


Maximum JV vt 


(C) 


140 


140 


170 


150 




(Cu) 


130 


120 


120 


130 




(Sn) 


110 


110 


120 


110 


LAI in accepted spill part 


(C) 


26231 


108215 


161331 


217899 




(Cu) 


27287 


87974 


175770 


90752 




(Sn) 


30029 


98078 


199209 


100872 


Fraction of triggers used 


(C) 


46 % 


42 % 


48 % 


43 % 




(Cu) 


46 % 


37 % 


32 % 


40 % 




(Sn) 


36 % 


32 % 


33 % 


34 % 


Accepted momentum 


(C) 


32483 


154984 


258338 


304993 


determination 


(Cu) 


37681 


156847 


374701 


209043 




(Sn) 


42949 


188994 


481436 


274700 


In kinematic region and 


(C) 


20508 


95999 


150444 


173077 


originating from target 


(Cu) 


23896 


99652 


229002 


122273 




(Sn) 


29090 


125864 


305214 


167137 


Negative particles 


(C) 


2873 


20328 


38892 


48699 




(Cu) 


3016 


18242 


52447 


31239 




(Sn) 


3352 


20721 


63846 


39395 


Positive particles 


(C) 


17635 


75671 


111552 


124378 




(Cu) 


20880 


81410 


176555 


91034 




(Sn) 


25738 


105143 


241368 


127742 


7r selected with PID 


(C) 


2661 


18513 


35115 


42994 




(Cu) 


2728 


16451 


46820 


27636 




(Sn) 


3100 


18762 


56697 


34595 


7r + selected with PID 


(C) 


5554 


28446 


47165 


54481 




(Cu) 


4403 


22087 


57218 


32234 




(Sn) 


4439 


23402 


64410 


38000 



and 12 GeV/c beams in p-C data, respectively. 

The double-differential cross-section for the production of a particle of type a can be expressed in the 
laboratory system as: 

d 2 CT„ 1 A „ . , , .... i ! ; 



d Pi d0j N pot N AP t ., V , vat Ja 1 3 ' 



where d ^ is expressed in bins of true momentum (pi), angle (8j) and particle type (a). 

The factor N ^ pt is the inverse of the number of target nuclei per unit area (A is the atomic mass, Na 
is the Avogadro number, p and t are the target density and thickness J3 The result is normalized to the 
number of incident protons on target N pot . 



2 We do not make a correction for the attenuation of the proton beam in the target, so that strictly speaking the 
cross-sections are valid for a Ai = 5% target. 



The 'raw yield' N-fj, is the number of particles of observed type a' in bins of reconstructed momentum 
(j>i>) and angle (9j>). These particles must satisfy the event, track and PID selection criteria. Although, 
owing to the stringent PID selection, the background from misidentified protons in the pion sample is 
small, the pion and proton raw yields (N"j,, for a' = 7r~,7r + ,p) have been measured simultaneously. 
This makes it possible to correct for the small remaining proton background in the pion data without 
prior assumptions concerning the proton production cross-section. 

The matrix Mj7 .,.,„,/ corrects for the efficiency and resolution of the detector. It unfolds the true 
variables ija from the reconstructed variables i j a with a Bayesian technique [22] and corrects the 
observed number of particles to take into account effects such as trigger efficiency, reconstruction efficiency, 
acceptance, absorption, pion decay, tertiary production, PID efficiency, PID misidentification and electron 
background. The method used to correct for the various effects is described in more detail in Ref. [T2] . 

In order to predict the population of the migration matrix element Mij a i'ji a i , the resolution, efficiency 
and acceptance of the detector are obtained from the Monte Carlo. This is accurate provided the Monte 
Carlo simulation describes these quantities correctly. Where some deviations from the control samples 
measured from the data are found, the data are used to introduce (small) ad hoc corrections to the Monte 
Carlo. 

Using the unfolding approach, possible known biases in the measurements are taken into account auto- 
matically as long as they are described by the Monte Carlo. For example the energy-loss of particles 
inside the target and material around the inner field cage is expressed as an average shift of the measured 
momentum distribution compared to the physical momentum. Known biases are therefore treated in the 
same way as resolution effects. In the experiment simulation, which is based on the GEANT4 toolkit [2"5] . 
the materials in the beam-line and the detector are accurately described as well as the relevant features 
of the detector response and the digitization process. The Monte Carlo simulation compares well with 
data, as shown in Ref. [12] . 

The absolute normalization of the result is calculated in first instance relative to the number of inci- 
dent beam particles accepted by the selection. After unfolding, the factor N ^ pt is applied. The beam 
normalization using down-scaled incident-proton triggers has uncertainties smaller than 2% for all beam 
momentum settings. 

The background due to interactions of the primary protons outside the target (called 'Empty target 
background') is measured using data taken without the target mounted in the target holder. Owing 
to the selection criteria which only accept events from the target region and the good definition of the 
interaction point this background is neglig iblc@ (< lCr 5 ). 

The effects of these uncertainties on the final results are estimated by repeating the analysis with the 
relevant input modified within the estimated uncertainty intervals. In many cases this procedure requires 
the construction of a set of different migration matrices. The correlations of the variations between the 
cross-section bins are evaluated and expressed in the covariance matrix. Each systematic error source is 
represented by its own covariance matrix. The sum of these matrices describes the total systematic error. 



4 Results 

The measured double-differential cross-sections for the production of ir + and 7r~ in the laboratory system 
as a function of the momentum and the polar angle for each incident beam momentum are shown in Fig. [2] 
and[3j respectively. The error bars shown are the square-roots of the diagonal elements in the covariance 
matrix, where statistical and systematic uncertainties are combined in quadrature. Correlations cannot 
be shown in the figures. The correlation of the statistical errors (introduced by the unfolding procedure) 
are typically smaller than 20% for adjacent momentum bins and smaller for adjacent angular bins. The 
correlations of the systematic errors are larger, typically 80% for adjacent bins. Tables with the results 

3 The background of interactions of the primary proton outside the target can be suppressed for large angle tracks 
measured in the TPC owing to the good resolution in z. This is contrary to the situation in the forward spectrometer where 
an interaction in the target cannot be distinguished from an interaction in upstream or downstream material |16l 117] . 
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of this analysis are also given in Appendix A. A discussion of the error evaluation is given below. The 
overall scale error (< 2%) is not shown. The measurements for the different beam momenta are overlaid in 
the same figure. For the 3 GeV/c data the point-to-point statistical error is larger than the systematic 
error, except for the lowest secondary momentum bin. Especially in the middle of the range (around 
400 MeV/c), the systematic error is small. Thus the fluctuations between the points are expected to 
be of statistical nature. In the first angular bins the momentum resolution is relatively large compared 
to the bin size such that the unfolding procedure tends to display statistical fluctuations over two bins. 
Since the treatment of the data sets taken with different beam momenta is identical, structures visible in 
the spectra at 3 GeV/c and not visible in the other data sets are not likely to be artefacts of the efficiency 
corrections. Overall trends in the shapes, i.e. structures extending over more than two bins are, however, 
to be considered significant. 

To better visualize the dependence on the incoming beam momentum, the same data averaged over the 
angular range (separately for the forward going and backward going tracks) covered by the analysis are 
shown separately for tt + and 7r~ in Fig. (?] The spectrum of pions produced in the backward direction is 
much steeper than that in the forward direction. 

The increase of the pion yield per proton is visible in addition to a change of spectrum towards higher 
momentum of the secondaries produced by higher momentum beams in the forward direction. 

The dependence of the integrated pion yields on the incident beam momentum is shown in Fig. [5] and 
compared with the p-Ta data taken with the same apparatus (Ref. |12j). The tt + and 7r~ yields integrated 
over the region 0.350 rad < 8 < 1.550 rad and 100 MeV/c < p < 700 MeV/c are shown in the left panel 
and the data integrated over the region 0.350 rad < 8 < 0.950 rad and 250 MeV/c < p < 500 McV/c in 
the right panel. The beam energy dependence of the yields is clearly different in the p-C data compared 
to the p-Ta data. The dependence in the p-C data is much more flat with a saturation of the yield 
between 8 GeV/c and 12 GeV/c (in both integration regions). The 7r + and ir~ production yields exhibit 
a different behaviour. The p-Cu and p-Sn data are more similar to the p-Ta data than the p-C data 
indicating a smooth transition between light and heavy target nuclei. 

The integrated ir~/n + ratio in the forward direction is displayed in Fig. [6] as a function of secondary 
momentum. The previously published p-Ta data are reproduced in addition to the measurements on the 
three target nuclei presented in this paper. In the covered part of the momentum range more 7r + 's are 
produced than 7r~'s. The ir~ /ir + ratio increases with increasing beam momentum and, depending on the 
beam momentum, a change of the sign of the slope of the ratio as a function of secondary momentum is 
visible in the p-C data. The latter feature is not present in the p-Cu, p-Sn and p-Ta [12] data. The ratio 
is closer to unity for the heavier target nuclei and a smaller variation with beam momentum is observed. 

The dependence of the integrated pion yields on the atomic number A is shown in Fig. [7] combining 
the results with the p-Ta data (Ref. [12] ) taken with the same apparatus and analysed using the same 
methods. The ir + yields integrated over the region 0.350 rad < 8 < 1.550 rad and 100 MeV/c < p < 
700 MeV/c are shown in the left panel and the ir~ data integrated over the same region in the right 
panel for four different beam momenta. One observes a smooth behaviour of the integrated yields. The 
A-dependence is slightly different for 7r~ and tt + production, the latter saturating earlier, especially at 
lower beam momenta. 

4.1 Systematic errors 

The uncertainties are reported in some detail in Table [2] for the carbon target data and summarized for 
the copper and tin target data in Table [31 One observes that only for the 3 GeV/c beam is the statistical 
error similar in magnitude to the systematic error, while the statistical error is negligible for the 8 GeV/c 
and 12 GeV/c beams. The statistical error is calculated by error propagation as part of the unfolding 
procedure. It takes into account that the unfolding matrix is obtained from the data themselvesQ and 
hence contributes also to the statistical error. This procedure almost doubles the statistical error, but 

4 The migration matrix is calculated without prior knowledge of the cross-sections, while the unfolding procedure deter- 
mined the unfolding matrix from the migration matrix and the distributions found in the data. 
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Figure 2: Double-differential cross-sections for ir + production in p-C, p-Cu and p-Sn interactions as 
a function of momentum displayed in different angular bins (shown in mrad in the panels). The results 
are given for four incident beam momenta (filled triangles: 3 GcV/c; open triangles: 5 GeV/c; filled 
rectangles: 8 GeV/c; open circles: 12 GcV/c). The error bars represent the combination of statistical 
and systematic uncertainties. 
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Figure 3: Double-differential cross-sections for ir~ production in p-C, p-Cu and p-Sn interactions as 
a function of momentum displayed in different angular bins (shown in mrad in the panels). The results 
are given for four incident beam momenta (filled triangles: 3 GcV/c; open triangles: 5 GeV/c; filled 
rectangles: 8 GeV/c; open circles: 12 GcV/c). The error bars represent the combination of statistical 
and systematic uncertainties. 
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Figure 4: Double-differential cross-sections for 7r + and 7r _ production in p-C, p-Cu and p-Sn interactions 
as a function of momentum averaged over the angular region covered by this experiment (shown in mrad) . 
The left panel of each pair shows forward production (350 mrad < 6 < 1550 mrad), while the right panel 
of each pair shows backward production (1550 mrad < 9 < 2150 mrad). The results are given for four 
incident beam momenta (filled triangles: 3 GeV/c; open triangles: 5 GeV/c; filled rectangles: 8 GeV/c; 
open circles: 12 GeV/c). The error bars obtained after summing the bins of the double-differential 
cross-sections take into account the correlations of the statistical and systematic uncertainties. 
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Figure 5: Left: The dependence on the beam momentum of the pion production yields in p-C, p-Cu, p- 
Sn, p-Ta interactions interactions integrated over the forward angular region (0.350 rad < 8 < 1.550 rad) 
and momentum (100 MeV/c < p < 700 MeV/c). Right: The dependence on the beam momentum of 
the pion production yields integrated over the region (0.350 rad < 9 < 0.950 rad and 250 MeV/c < p < 
500 MeV/c) with the same meaning of the symbols. The results are given in arbitrary units, with a 
consistent scale between the left and right panel. Although the units are indicated as "arbitrary" , for the 
largest region (left panel), the yield is expressed as d 2 cr/dpdil in mb/(GeV/c sr). For the smaller region 
(left panel) the same normalization is chosen, but now scaled with the relative bin size to show visually 
the correct ratio of number of pions produced in this kinematical region with respect to the yield in the 
larger kinematical region. Data points for different target nuclei and equal momenta are slightly shifted 
horizontally with respect to each other to increase the visibility. 



avoids an important systematic error which would otherwise be introduced by assuming a cross-section 
model a priori to calculate the corrections. 

The largest systematic error corresponds to the uncertainty in the absolute momentum scale, which was 
estimated to be around 3% using elastic scattering (sec detailed discussion in [12 ). It is difficult to better 
constrain the absolute momentum scale, since it depends on the knowledge of the beam momentum 
(known to 1%) and the measurement of the forward scattering angle in the elastic scattering interaction. 
At low momentum in the relatively small angle forward direction the uncertainty in the subtraction of the 
electron and positron background due to ir° production is dominant. This uncertainty is split between 
the variation in the shape of the tt° spectrum and the normalization using the recognized electrons. The 
assumption is made that the ir° spectrum is similar to the spectrum of charged pions. Initial ir~ and ir + 
spectra are obtained in an analysis without 7r° subtraction. The ir~ spectra are then used in the MC for 
the 7T° distributions. A full simulation of the production and decay into 7's with subsequent conversion 
in the detector materials is used to predict the background electron and positron tracks. In the region 
below 120 MeV/c a large fraction of the electrons can be unambiguously identified. These tracks are 
used as relative normalization between data and MC. The remaining background is then estimated from 
the distributions of the simulated electron and positron tracks which are accepted as pion tracks with 
the same criteria as used to select the data. These normalized distributions are subtracted from the 
data before the unfolding procedure is applied. Uncertainties in the assumption of the ir° spectrum are 
taken into account by an alternative assumption that their spectrum follows the average of the tt~ and 7r + 
distribution. An additional systematic error of 10% is assigned to the normalization of the ir° subtraction 
using the identified electrons and positrons. 



14 



HARP p-C 77-/tt + 



HARP p-Cu tt'/tt* 



1.5 



0.5 



-O 12 
" ■ 8 

A 5 



1.5 



0.5 



2 

350-1550 mrad 



1.5 




hAi^A-A; 



0.2 0.4 
P [GeV/c] 
HARP p-Sn ttV" + 



0.5 



350-1550 mrad 



■O 12 
■ 8 

A 5 




0.6 



0.2 0.4 
P [GeV/c] 
HARP p-Ta 7T/7T+ 



0.6 




0.2 0.4 

P [GeV/c] 



0.2 0.4 0.6 

P [GeV/c] 



Figure 6: The ratio of the differential cross-sections for 7r~ and n + production in p-C, p-Cu, p-Sn 
and p-Ta interactions as a function of secondary momentum integrated over the forward angular region 
(shown in mrad). The results are given for four incident beam momenta (filled triangles: 3 GeV/c; open 
triangles: 5 GeV/c; filled rectangles: 8 GeV/c; open circles: 12 GeV/c). 
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Figure 7: The dependence on the atomic number A of the pion production yields in p-C, p-Cu, p-Sn, 
p-Ta interactions interactions integrated over the forward angular region (0.350 rad < 9 < 1.550 rad) and 
momentum (100 MeV/c < p < 700 MeV/c). The results are given in arbitrary units, with a consistent 
scale between the left and right panel. The vertical scale used in this figure is consistent with the one in 
Fig. El 

The target region definition and the uncertainty in the PID efficiency and background from tertiaries are 
of similar size and are not negligible. Relatively small errors are introduced by the uncertainties in the 
absorption correction, absolute knowledge of the angular and the momentum resolution. The correction 
for tertiaries (particles produced in secondary interactions) is relatively large at low momenta and large 
angles. As expected, this region is most affected by this component. 

As already mentioned above, the overall normalization has an uncertainty of 2%, and is not reported 
in the table. It is mainly due to the uncertainty in the efficiency that beam protons counted in the 
normalization actually hit the target, with smaller components from the target density and beam particle 
counting procedure. 



4.2 Comparisons with earlier data 

Very few pion production data sets are available in the literature for p-C, p-Cu and p-Sn interactions in 
this energy region. Our data can be compared with results from Ref. [24] and [25] where measurements 
of 7r~ production are reported in 4.2 GeV/c and 10 GeV/c p-C interactions, respectively. The total 
number of 7r~ observed in the above references is about 1300 (5650) in the 4.2(10) GeV/c data. In the 
papers cited above no tables of the double differential cross-sections were provided, the measurements 
being given in parametrized and graphical form only. The authors of Ref. [24] and [25] give the results as 
a simple exponential in the invariant cross-section: ^ ^Lf , where E and p are the energy and momentum 
of the produced particle, respectively, and A the atomic number of the target nucleu£0. Unfortunately, no 
absolute normalization is given numerically. To provide a comparison with these data, the parametrization 

5 their spectra are parametrized in each angular bin with a function of the form / _ = c exp (— T/To), where T is the 
kinetic energy of the produced particle and To is given by To = T'/(l — /3 cos 8). For the 4.2 GeV/c data the values of the 
parameters are X" = (0.089 ± 0.006) GeV/c and j3 = 0.77 ± 0.04 and T" = (0.100 ± 0.002) GeV/c and j3 = 0.81 ± 0.02 for 
the 10 GeV/c data. 
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Table 2: Contributions to the experimental uncertainties for the thin carbon target data. The numbers 
represent the uncertainty in percent of the cross-section integrated over the angle and momentum region 
indicated. The overall normalization has an uncertainty of 2%, and is not reported in the table. 
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Table 3: Summary of experimental uncertainties for the copper (tin) analysis. The numbers represent 
the uncertainty in percent of the cross-section integrated over the angle and momentum region indicated. 
The overall normalization has an uncertainty of 2%, and is not reported in the table. 
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was integrated over the angular bins used in our analysis and with an arbitrary overall normalization 
overlaid to our results. We compare the 4.2 GeV/c parametrization of Ref. [53] with our 5 GeV/c 
data and the Ref. [55] parametrization with our 12 GeV/c data. In the comparison with the 4.2 GeV/c 
parametrization the normalization c is a simple constant, while for the 10 GeV/c parametrization a smooth 
^-dependence consistent with a graphical analysis of Ref. [25 was used. Thus only the comparison of 
the slopes with secondary momentum can be considered significant. Since the 8 GeV/c and 12 GeV/c 
p-C results are very similar, the lack of data with an exactly equal beam momentum does not play an 
important role. The results of this comparison are shown in Fig. [8] The shaded band gives the excursion 
of the parametrization due to the error in the slope parameters (±2c) with an additional assumed 10% 
error on the absolute scale. The latter additional error takes into account the fact that the errors on the 
slopes fitted to the individual angular bins in the cited data are at least a factor of two larger than in 
the exponential slope obtained from their global parametrization. The agreement of our data with the 
simple parametrization is good. To judge the quality of the comparison, one should keep in mind that 
the statistics of Ref. [53] and [55] is much smaller (1300 7r~ and 5650 7r~, respectively) than the statistics 
of the 7r~ samples in our 5 GeV/c and 12 GeV/c data (18,000 and 43,000 7r~, respectively). The errors 
on the slopes fitted to the individual angular bins in the cited data are at least a factor of two larger than 
in the exponential slope obtained from their global parametrization. The bands in the figure extend over 
the region where data from Ref. [24] and Ref. [25] are available. 

Our p-C and p-Cu data can also be compared with 7r + and 7r~ production measurements taken with 
12 GeV/c incident protons from Ref. [26]. These data were taken with a magnetic spectrometer and only 
measurements at 90 degrees from the initial proton direction are available. The statistical point-to-point 
errors are quoted to be 3%, while the overall normalization has a 30% uncertainty due to the knowledge 
of the acceptance. In Fig. [5] their p-C data are shown together with the p-C data reported in this paper. 
The filled boxes show the data directly from Ref. [55], while the open boxes are scaled with a factor 0.72. 
This factor was defined by scaling the average of the tt~ and 7r + data from Ref. [26, at 179 MeV/c and 
242 MeV/c to the HARP data averaged over the same region. The scale factor is within one standard 
deviation of the systematic normalization uncertainty of Ref. [26j . The latter data set compares well with 
the data described in this paper (filled circles) in the angular region 1.35 rad < 9 <1.55 rad at the same 
proton beam momentum. 

In Fig. [TU] the p-Cu data of Ref. [55] are shown together with the p-Cu results reported in this paper. 
The filled boxes show the data directly from Ref. [55J, while the open boxes are scaled with a factor 
0.91. This factor was defined in a similar procedure as described for the p-C data. The scale factor is 
very close to unity and well within one standard deviation of the systematic normalization uncertainty 
of Ref. [26] . Like for the p-C data, the full error bars including the 30% scale uncertainty of Ref. [26] are 
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Figure 8: Comparison of the HARP p-C results with data from Ref. [24] and [25) . The left panel shows 
the comparison of the parametrization of the 4.2 GeV/c data of Ref. [24j with the 5 GeV/c data reported 
here; the right panel shows the comparison of the 10 GeV/c parametrization of [25 with the 12 GeV/c 
data. The absolute normalization of the parametrization was fixed to the data in both cases. The 
band shows the range allowed by varying the slope parameters given by [24] and [25] with two standard 
deviation and a 10% variation on the absolute scale. The angular ranges are shown in mrad in the panels. 



drawn, for the scaled data only their quoted statistical error of 3% is shown. The agreement of the two 
data sets is excellent. The fact that the two scale factors are different may be due to the fact that the 
scale uncertainty in Ref. [26] holds separately for data sets taken with different target nuclei. 

Available data at 12.3 GeV/c from the E910 experiment [27] are in reasonable agreement with our p-Cu 
results as shown in Fig. 1111 In order to take into account the different angular binnings which prevent 
a direct comparison, a Sanford-Wang parametrization is fitted to our data. The fit is performed to the 
data redefined as d 2 a^ /dpdfl . An area between two parametrizations is defined which contains our data 
points as shown in Fig. [Tl] (top panels). It is visible that the parametrization is not a perfect description 
to our data. Therefore, we define a band of ±15% around the best fit which contains almost all the 
HARP data points. The same parametrizations are then displayed in the binning of E910. While the 
shape of the distributions are similar for both ir + and ir~ in HARP and E910 data sets, the absolute 
normalizations disagree by 5%-10%. For the individual data sets the systematic errors are between 5% 
and 10% depending on the range of secondary momentum. Since these errors are correlated between 
bins, the discrepancy in the tt + and n~ data separately are of the order of one standard deviation. 
However, the effects are opposite in tt + and 7r~, giving a 15% difference in the tt + /it~ ratio between two 
experiments which is of the order of two standard deviations. This effect may point to an underestimation 
of systematic effects on the absolute normalization in one of the experiments or in the PID efficiency. 
Part of the difference is also due to an imperfect parametrization of our data sample. Owing to the 
symmetry of the HARP TPC, including its trigger counter, we do not expect a large systematic error in 
the HARP data between ir + and tt~ production cross-sections. 
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Figure 9: Comparison of the HARP results with ir + and ir~ production data at 90 degrees from Ref. [2"6] 
taken with 12 GeV/c protons. The left panel shows the comparison of the tt + production data of Ref. [26] 
with the data reported here; the right panel shows the comparison with the 7r~ production data. The 
smaller filled boxes show the data directly from Ref. [26], while the open boxes are scaled as explained 
in the text. The latter data set compares well with the data described in this paper (filled circles) in the 
angular region 1.35 rad < 8 <1.55 rad. 



5 Summary and Conclusions 



An analysis of the production of pions at large angles with respect to the beam direction for protons of 
3 GeV/c, 5 GeV/c, 8 GeV/c and 12 GeV/c impinging on thin (5% interaction length) carbon, copper 
and tin targets is described. The secondary pion yield is measured in a large angular and momentum 
range and double-differential cross-sections are obtained. A detailed error estimation has been discussed. 
Results on the dependence of pion production on the target atomic number A are also presented. 

The use of a single detector for a range of beam momenta makes it possible to measure the dependence 
of the pion yield on the secondary particle momentum and emission angle 9 with high precision. The A 
dependence of the cross-section can be studied using the combination of the present data with the data 
obtained with tantalum [12j . 

Very few pion production measurements in this energy range are reported in the literature. The only 
comparable results found in the literature agrees with the analysis described in this paper. Hadronic 
production models describing this energy range can now be compared with our new results and, if needed, 
improved. Data taken with different target materials and beam momenta will be presented in subsequent 
papers. 
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Figure 11: Comparison of the HARP data with tt + and 7r _ production data from Ref. [57] taken with 
12.3 GeV/c protons. The top panels show a parametrization of the 7r + (left) and ir~ (right) production 
data described in this paper. The data have been normalized to represent d 2 <r v /dpdfl . The shaded band 
represents the area between two parametrization which contain the data points. The bottom panels show 
the comparison of the same parametrization, now binned according to the E910 data. The bottom left 
(right) panel shows the ir + production data of Ref. [27j . The angular regions are indicated in mrad 
in the upper right-hand corner of each plot. 
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A Cross-section data 



Table 4: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 ^ /{dpdO) for carbon. Each row refers to a different (p m i n < p < Pmax,#min — ^ < #max) bin, 
where p and 8 are the pion momentum and polar angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 
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0.114±0.015 






0.15 


0.20 


0.114±0.016 


0.131±0.009 


0.141±0.009 


0.131±0.008 






0.20 


0.25 


0.093±0.014 


0.088±0.008 


0.088±0.008 


0.094±0.006 






0.25 


0.30 


0.035±0.010 


0.045±0.005 


0.052±0.004 


0.053±0.006 






0.30 


0.35 


0.013±0.005 


0.026±0.004 


0.034±0.004 


0.039±0.003 






0.35 


0.40 


0.005±0.003 


0.014±0.003 


0.018±0.003 


0.025±0.004 






0.40 


0.45 


0.003±0.002 


0.006±0.002 


0.010±0.002 


0.012±0.003 






0.45 


0.50 


0.002±0.002 


0.003±0.001 


0.006±0.001 


0.006±0.002 


1.95 


2.15 


0.10 


0.15 


0.073±0.015 


0.087±0.012 


0.110±0.013 


0.095±0.012 






0.15 


0.20 


0.096±0.016 


0.094±0.008 


0.106±0.008 


0.108±0.008 






0.20 


0.25 


0.059±0.012 


0.059±0.006 


0.057±0.005 


0.063±0.006 






0.25 


0.30 


0.026±0.009 


0.027±0.004 


0.036±0.004 


0.032±0.004 






0.30 


0.35 


0.012±0.005 


0.017±0.002 


0.023±0.002 


0.016±0.002 






0.35 


0.40 


0.004±0.003 


0.014±0.002 


0.014±0.002 


0.011±0.002 






0.40 


0.45 


0.003±0.003 


0.005±0.002 


0.006±0.002 


0.007±0.001 






0.45 


0.50 


0.002±0.002 


0.002±0.001 


0.002±0.001 


0.004±0.001 
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Table 5: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d 2 a n /(dpd9) for carbon. Each row refers to a different (p m j n < p < Pmax, #min < # < #max) bin, 
where p and 9 are the pion momentum and polar angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 





$max 




Pmax 




d 2 a^ 


/(dpd9) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c rad)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.35 


0.55 


0.10 


0.15 


0.015±0.020 


08±0 04 

V 7 . V / v. 7 — 1 — V 7 . V / 1 


13±0 06 

V 7 . 1 1 1 — l — V 7 ■ V / \7 


0.10±0.05 






0.15 


0.20 


0.036±0.023 


076±0 025 


0.113±0.025 


0.115±0.028 






0.20 


0.25 


0.053±0.020 


0.109±0.018 


0.162±0.023 


0.172±0.029 






0.25 


0.30 


0.066±0.015 


0.130±0.016 


0.179±0.017 


0.220±0.020 






0.30 


0.35 


0.066±0.014 


0.133±0.013 


0.204±0.019 


0.229±0.016 






0.35 


0.40 


0.069±0.014 


0.146±0.013 


0.193±0.013 


0.216±0.012 






0.40 


0.45 


0.082±0.012 


0.126±0.009 


0.193±0.014 


0.212±0.014 






0.45 


0.50 


0.078±0.012 


0.113±0.007 


0.177±0.011 


0.223±0.012 






0.50 


0.60 


0.045±0.008 


0.120±0.008 


0.186±0.012 


0.222±0.012 






0.60 


0.70 


0.050±0.009 


0.118±0.010 


0.176±0.015 


0.212±0.017 






0.70 


0.80 


0.052±0.010 


0.100±0.011 


0.152±0.017 


0.194±0.020 


0.55 


0.75 


0.10 


0.15 


0.040±0.026 


0.049±0.022 


0.053±0.024 


0.090±0.030 






0.15 


0.20 


0.070±0.018 


0.095±0.019 


0.140±0.019 


0.143±0.017 






0.20 


0.25 


0.089±0.017 


136±0 014 

V 7 . 1 1 1 V7 — 1 — V 7 . V / J 1 


166±0 017 

V / ■ 1 V7 V7 — 1 — V 7 . V / 1 1 


0.177±0.015 






0.25 


0.30 


0.072±0.013 


0.147±0.012 


0.187±0.013 


0.211±0.017 






0.30 


0.35 


0.098±0.016 


0.128±0.009 


0.166±0.012 


0.206±0.012 






0.35 


0.40 


0.081±0.013 


0.131±0.011 


0.170±0.010 


0.184±0.009 






0.40 


0.45 


0.077±0.011 


0.111±0.009 


0.165±0.009 


0.181±0.009 






0.45 


0.50 


0.072±0.010 


0.089±0.006 


0.158±0.009 


0.178±0.009 






0.50 


0.60 


0.049±0.008 


0.097±0.006 


0.145±0.009 


0.164±0.009 






0.60 


0.70 


0.035±0.007 


0.081±0.008 


0.119±0.010 


0.133±0.012 






0.70 


0.80 


0.029±0.007 


0.057±0.010 


0.099±0.013 


0.104±0.014 


0.75 


0.95 


0.10 


0.15 


0.048±0.018 


063±0 016 

V 7 . V / \7 1 7 — 1 — V 7 . V / 1 \7 


076±0 016 

V7 . V 7 i V7 — 1 — V 7 . V / 1 V 7 


0.087±0.017 






0.15 


0.20 


0.064±0.015 


0.124±0.014 


0.175±0.018 


0.187±0.019 






0.20 


0.25 


0.065±0.014 


0.125±0.011 


0.184±0.015 


0.187±0.012 






0.25 


0.30 


0.087±0.014 


0.124±0.011 


0.177±0.011 


0.180±0.013 






0.30 


0.35 


0.062±0.010 


0.119±0.009 


0.148±0.009 


0.173±0.009 






0.35 


0.40 


0.065±0.011 


0.110±0.008 


0.141±0.010 


0.142±0.007 






0.40 


0.45 


0.049±0.007 


0.093±0.007 


0.128±0.007 


0.132±0.006 






0.45 


0.50 


0.049±0.008 


0.070±0.006 


0.107±0.007 


0.119±0.005 






0.50 


0.60 


0.038±0.007 


0.066±0.004 


0.087±0.005 


0.101±0.006 






0.60 


0.70 


0.018±0.006 


0.055±0.006 


0.074±0.007 


0.074±0.007 


0.95 


1.15 


0.10 


0.15 


0.038±0.013 


0.068±0.012 


0.094±0.016 


0.101±0.013 






0.15 


0.20 


0.071±0.014 


0.129±0.011 


0.141±0.010 


0.138±0.014 






0.20 


0.25 


0.074±0.013 


0.130±0.011 


0.137±0.012 


0.167±0.012 






0.25 


0.30 


0.094±0.014 


0.116±0.008 


0.134±0.009 


0.152±0.009 






0.30 


0.35 


0.057±0.011 


0.093±0.006 


0.126±0.010 


0.129±0.007 






0.35 


0.40 


0.042±0.007 


0.079±0.005 


0.095±0.006 


0.099±0.006 






0.40 


0.45 


0.035±0.007 


0.065±0.004 


0.078±0.004 


0.084±0.004 






0.45 


0.50 


0.026±0.006 


0.053±0.004 


0.069±0.004 


0.077±0.004 






0.50 


0.60 


0.014±0.004 


0.039±0.004 


0.051±0.004 


0.061±0.004 
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^miri 


$max 


-Pmin 


Pmax 




d 2 a^ 


/(dpd8) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c rad)) 












3 GeV/c 


5 GcV/c 


8 GcV/c 


12 GeV/c 


1.15 


1.35 


0.10 


0.15 


0.061±0.015 


0.070±0.010 


0.103±0.012 


0.106±0.013 






0.15 


0.20 


0.124±0.019 


0.125±0.013 


0.139±0.010 


0.158±0.011 






0.20 


0.25 


0.070±0.012 


0.121±0.009 


0.135±0.010 


0.135±0.008 






0.25 


0.30 


0.056±0.011 


0.089±0.007 


0.111±0.008 


0.119±0.007 






0.30 


0.35 


0.036±0.010 


0.073±0.006 


0.082±0.005 


0.099±0.005 






0.35 


0.40 


0.016±0.004 


0.057±0.005 


0.068±0.004 


0.078±0.004 






0.40 


0.45 


0.013±0.004 


0.037±0.004 


0.056±0.004 


0.059±0.004 






0.45 


0.50 


0.009±0.003 


0.027±0.003 


0.043±0.004 


0.050±0.004 


1.35 


1.55 


0.10 


0.15 


0.046±0.013 


0.071±0.009 


0.095±0.011 


0.108±0.013 






0.15 


0.20 


0.093±0.016 


0.105±0.010 


0.129±0.012 


0.134±0.010 






0.20 


0.25 


0.051±0.011 


0.093±0.008 


0.110±0.007 


0.116±0.007 






0.25 


0.30 


0.052±0.010 


0.078±0.007 


0.087±0.006 


0.079±0.005 






0.30 


0.35 


0.033±0.007 


0.045±0.005 


0.054±0.005 


0.065±0.004 






0.35 


0.40 


0.026±0.006 


0.032±0.003 


0.046±0.003 


0.048±0.004 






0.40 


0.45 


0.024±0.006 


0.024±0.003 


0.037±0.003 


0.032±0.003 






0.45 


0.50 


0.015±0.005 


0.018±0.002 


0.026±0.003 


0.026±0.002 


1.55 


1.75 


0.10 


0.15 


0.077±0.016 


0.075±0.012 


0.095±0.012 


0.090±0.011 






0.15 


0.20 


0.094±0.016 


0.117±0.010 


0.117±0.008 


0.127±0.009 






0.20 


0.25 


0.045±0.010 


0.065±0.006 


0.077±0.005 


0.087±0.007 






0.25 


0.30 


0.022±0.007 


0.057±0.006 


0.061±0.005 


0.061±0.005 






0.30 


0.35 


0.010±0.004 


0.036±0.004 


0.048±0.004 


0.043±0.004 






0.35 


0.40 


0.006±0.003 


0.028±0.003 


0.038±0.003 


0.032±0.003 






0.40 


0.45 


0.008±0.004 


0.018±0.003 


0.026±0.003 


0.024±0.002 






0.45 


0.50 


0.007±0.004 


0.012±0.002 


0.017±0.003 


0.017±0.002 


1.75 


1.95 


0.10 


0.15 


0.077±0.016 


0.074±0.010 


0.087±0.011 


0.069±0.008 






0.15 


0.20 


0.082±0.015 


0.095±0.008 


0.101±0.007 


0.103±0.008 






0.20 


0.25 


0.025±0.007 


0.067±0.007 


0.069±0.006 


0.070±0.005 






0.25 


0.30 


0.021±0.007 


0.040±0.005 


0.046±0.005 


0.051±0.004 






0.30 


0.35 


0.013±0.005 


0.024±0.003 


0.026±0.003 


0.032±0.003 






0.35 


0.40 


0.011±0.006 


0.016±0.002 


0.020±0.002 


0.022±0.002 






0.40 


0.45 


0.003±0.003 


0.011±0.002 


0.018±0.002 


0.017±0.002 






0.45 


0.50 


0.001±0.001 


0.007±0.002 


0.013±0.002 


0.011±0.002 


1.95 


2.15 


0.10 


0.15 


0.040±0.010 


0.058±0.008 


0.066±0.007 


0.069±0.007 






0.15 


0.20 


0.071±0.014 


0.081±0.009 


0.081±0.008 


0.080±0.007 






0.20 


0.25 


0.048±0.012 


0.044±0.005 


0.064±0.005 


0.062±0.005 






0.25 


0.30 


0.012±0.006 


0.025±0.004 


0.038±0.004 


0.038±0.003 






0.30 


0.35 


0.003±0.002 


0.010±0.002 


0.025±0.003 


0.022±0.003 






0.35 


0.40 


0.002±0.002 


0.007±0.001 


0.014±0.002 


0.015±0.002 






0.40 


0.45 


0.002±0.002 


0.006±0.001 


0.008±0.002 


0.012±0.002 






0.45 


0.50 


0.003±0.003 


0.004±0.001 


0.005±0.001 


0.007±0.001 
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Table 6: HARP results for the double-differential ir + production cross-section in the laboratory system, 
cPa n /(dpdO) for copper. Each row refers to a different (p m i n < p < Pmax,^min < " < "max) bin, 
where p and 9 are the pion momentum and polar angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 





"max 


.Pmin 


Pmax 






/{dpdff) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c rad)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.35 


0.55 


0.10 


0.15 


0.31±0.14 


0.41±0.17 


46±0 20 


0.69±0.29 






0.15 


0.20 


0.17±0.05 


41±0 09 

v / ■ i i — i — v / ■ \ > • ' 


63±0 11 


0.84±0.14 






0.20 


0.25 


0.22±0.07 


0.61±0.06 


0.82±0.07 


0.83±0.07 






0.25 


0.30 


0.52±0.08 


0.75±0.07 


0.96±0.07 


1.13±0.10 






0.30 


0.35 


0.43±0.06 


0.79±0.06 


1.08±0.09 


1.22±0.10 






0.35 


0.40 


0.42±0.05 


0.75±0.05 


1.13±0.06 


1.41±0.09 






0.40 


0.45 


0.33±0.04 


0.77±0.05 


1.09±0.06 


1.19±0.08 






0.45 


0.50 


0.41±0.06 


0.70±0.05 


1.07±0.06 


1.17±0.09 






0.50 


0.60 


0.44±0.05 


0.69±0.04 


1.08±0.06 


1.16±0.08 






0.60 


0.70 


0.25±0.05 


0.56±0.06 


0.92±0.09 


1.18±0.11 






0.70 


0.80 


0.13±0.03 


0.37±0.07 


0.72±0.10 


0.92±0.14 


0.55 


0.75 


0.10 


0.15 


0.34±0.11 


0.37±0.11 


0.48±0.14 


0.53±0.15 






0.15 


0.20 


0.48±0.07 


0.64±0.06 


0.82±0.08 


0.80±0.08 






0.20 


0.25 


0.44±0.06 


76±0 06 

V / ■ 1 \ } — 1 — V / ■ \ J \ / 


1.04±0.09 


1.18±0.11 






0.25 


0.30 


0.50±0.06 


0.82±0.08 


1.01±0.06 


1.29±0.09 






0.30 


0.35 


0.51±0.06 


0.89±0.06 


1.06±0.08 


1.12±0.07 






0.35 


0.40 


0.47±0.05 


0.75±0.04 


1.06±0.06 


1.19±0.08 






0.40 


0.45 


0.38±0.04 


0.64±0.04 


0.97±0.05 


1.12±0.06 






0.45 


0.50 


0.33±0.04 


0.57±0.03 


0.91±0.04 


0.99±0.06 






0.50 


0.60 


0.26±0.03 


0.47±0.04 


0.75±0.05 


0.85±0.05 






0.60 


0.70 


0.15±0.03 


0.33±0.04 


0.49±0.06 


0.63±0.07 






0.70 


0.80 


0.06±0.02 


0.20±0.04 


0.30±0.06 


0.43±0.08 


0.75 


0.95 


0.10 


0.15 


0.41±0.10 


43±0 10 

v / ■ i 'j — i — v / ■ ivy 


0.52±0.12 


0.63±0.13 






0.15 


0.20 


0.54±0.06 


0.73±0.06 


0.89±0.06 


1.04±0.10 






0.20 


0.25 


0.56±0.07 


0.81±0.06 


1.06±0.09 


1.23±0.08 






0.25 


0.30 


0.53±0.06 


0.73±0.05 


1.03±0.06 


1.06±0.07 






0.30 


0.35 


0.38±0.05 


0.66±0.04 


0.92±0.06 


1.04±0.06 






0.35 


0.40 


0.27±0.03 


0.54±0.03 


0.81±0.04 


0.92±0.05 






0.40 


0.45 


0.24±0.03 


0.45±0.03 


0.63±0.03 


0.82±0.04 






0.45 


0.50 


0.21±0.03 


0.37±0.02 


0.55±0.03 


0.67±0.04 






0.50 


0.60 


0.10±0.02 


0.27±0.02 


0.43±0.03 


0.51±0.04 






0.60 


0.70 


0.05±0.01 


0.16±0.03 


0.26±0.04 


0.32±0.05 


0.95 


1.15 


0.10 


0.15 


0.35±0.07 


0.62±0.11 


0.63±0.11 


0.77±0.15 






0.15 


0.20 


0.43±0.06 


0.81±0.05 


1.02±0.07 


1.06±0.08 






0.20 


0.25 


0.41±0.05 


0.59±0.05 


0.97±0.05 


1.03±0.07 






0.25 


0.30 


0.37±0.05 


0.63±0.05 


0.81±0.05 


0.82±0.05 






0.30 


0.35 


0.28±0.04 


0.51±0.04 


0.64±0.04 


0.72±0.04 






0.35 


0.40 


0.17±0.02 


0.35±0.02 


0.52±0.03 


0.59±0.03 






0.40 


0.45 


0.16±0.03 


0.29±0.02 


0.39±0.03 


0.47±0.03 






0.45 


0.50 


0.11±0.02 


0.23±0.02 


0.29±0.02 


0.39±0.03 






0.50 


0.60 


0.05±0.01 


0.15±0.02 


0.19±0.02 


0.25±0.03 
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^min 
(rad) 


(rad) 


-Pmin 
(GeV/c) 


Pmax 
(GeV/c) 




d 2 a* + /{dpd6) 
(barn/(GeV/c rad)) 












3 GeV/c 


5 GeV/c 


8 GcV/c 


12 GeV/c 


1.15 


1.35 


0.10 


0.15 


0.43±0.09 


0.60±0.13 


0.66±0.14 


0.67±0.14 






0.15 


0.20 


0.51±0.06 


0.80±0.05 


0.96±0.07 


0.96±0.10 






0.20 


0.25 


0.42±0.05 


0.56±0.04 


0.87±0.05 


1.04±0.06 






0.25 


0.30 


0.37±0.05 


0.44±0.03 


0.60±0.04 


0.74±0.04 






0.30 


0.35 


0.18±0.03 


0.31±0.03 


0.42±0.03 


0.51±0.03 






0.35 


0.40 


0.11±0.02 


0.23±0.02 


0.33±0.02 


0.37±0.03 






0.40 


0.45 


0.08±0.02 


0.16±0.02 


0.25±0.02 


0.27±0.02 






0.45 


0.50 


0.05±0.01 


0.11±0.01 


0.18±0.02 


0.19±0.02 


1.35 


1.55 


0.10 


0.15 


0.43±0.12 


0.64±0.15 


0.69±0.17 


0.75±0.16 






0.15 


0.20 


0.46±0.06 


0.73±0.06 


0.89±0.08 


0.88±0.09 






0.20 


0.25 


0.30±0.04 


0.48±0.04 


0.70±0.04 


0.79±0.05 






0.25 


0.30 


0.18±0.03 


0.34±0.03 


0.45±0.04 


0.50±0.04 






0.30 


0.35 


0.13±0.02 


0.24±0.02 


0.30±0.02 


0.36±0.03 






0.35 


0.40 


0.09±0.02 


0.15±0.01 


0.22±0.02 


0.28±0.02 






0.40 


0.45 


0.06±0.01 


0.09±0.01 


0.14±0.01 


0.20±0.02 






0.45 


0.50 


0.03±0.01 


0.06±0.01 


0.10±0.01 


0.12±0.02 


1.55 


1.75 


0.10 


0.15 


0.52±0.13 


0.63±0.16 


0.79±0.18 


0.75±0.18 






0.15 


0.20 


0.39±0.05 


0.77±0.06 


0.79±0.06 


0.80±0.07 






0.20 


0.25 


0.30±0.04 


0.40±0.04 


0.60±0.04 


0.62±0.05 






0.25 


0.30 


0.25±0.04 


0.23±0.02 


0.32±0.03 


0.34±0.03 






0.30 


0.35 


0.11±0.03 


0.16±0.02 


0.23±0.02 


0.23±0.02 






0.35 


0.40 


0.04±0.01 


0.10±0.01 


0.17±0.02 


0.18±0.02 






0.40 


0.45 


0.03±0.01 


0.07±0.01 


0.10±0.01 


0.12±0.02 






0.45 


0.50 


0.02±0.01 


0.05±0.01 


0.06±0.01 


0.07±0.01 


1.75 


1.95 


0.10 


0.15 


0.59±0.11 


0.54±0.10 


0.66±0.12 


0.66±0.12 






0.15 


0.20 


0.46±0.05 


0.61±0.04 


0.67±0.04 


0.68±0.05 






0.20 


0.25 


0.36±0.05 


0.34±0.03 


0.41±0.03 


0.43±0.04 






0.25 


0.30 


0.08±0.03 


0.16±0.02 


0.21±0.02 


0.25±0.02 






0.30 


0.35 


0.03±0.01 


0.11±0.01 


0.12±0.01 


0.14±0.02 






0.35 


0.40 


0.02±0.01 


0.08±0.01 


0.09±0.01 


0.09±0.01 






0.40 


0.45 


0.02±0.01 


0.05±0.01 


0.05±0.01 


0.05±0.01 






0.45 


0.50 


0.01±0.01 


0.02±0.01 


0.03±0.01 


0.03±0.01 


1.95 


2.15 


0.10 


0.15 


0.31±0.07 


0.45±0.07 


0.49±0.07 


0.47±0.08 






0.15 


0.20 


0.41±0.05 


0.44±0.04 


0.52±0.03 


0.50±0.04 






0.20 


0.25 


0.16±0.04 


0.22±0.02 


0.28±0.02 


0.28±0.02 






0.25 


0.30 


0.06±0.02 


0.10±0.02 


0.14±0.02 


0.20±0.02 






0.30 


0.35 


0.03±0.01 


0.06±0.01 


0.09±0.01 


0.11±0.02 






0.35 


0.40 


0.02±0.01 


0.05±0.01 


0.04±0.01 


0.06±0.01 






0.40 


0.45 


0.01±0.01 


0.02±0.01 


0.03±0.01 


0.04±0.01 






0.45 


0.50 






0.02±0.01 


0.02±0.01 
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Table 7: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa n /(dpd9) for copper. Each row refers to a different (p m j n < p < Pmax,"min < " < "max) bin, 
where p and 9 are the pion momentum and polar angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 





"max 


.Pmin 


Pmax 




d 2 a^ 


/(dpd9) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c rad)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.35 


0.55 


0.10 


0.15 


0.19±0.11 


45±0 19 


54±0 23 


0.55±0.27 






0.15 


0.20 


0.12±0.07 


44±0 09 


0.61±0.11 


0.81±0.15 






0.20 


0.25 


0.19±0.05 


0.48±0.06 


0.76±0.08 


0.96±0.09 






0.25 


0.30 


0.20±0.06 


0.49±0.05 


0.80±0.06 


0.94±0.07 






0.30 


0.35 


0.25±0.05 


0.45±0.04 


0.78±0.05 


0.92±0.07 






0.35 


0.40 


0.16±0.03 


0.43±0.04 


0.76±0.04 


0.89±0.06 






0.40 


0.45 


0.13±0.02 


0.40±0.03 


0.72±0.04 


0.79±0.04 






0.45 


0.50 


0.21±0.05 


0.41±0.03 


0.70±0.04 


0.79±0.05 






0.50 


0.60 


0.17±0.03 


0.41±0.03 


0.70±0.04 


0.79±0.05 






0.60 


0.70 


0.11±0.02 


0.32±0.03 


0.64±0.05 


0.69±0.06 






0.70 


0.80 


0.10±0.02 


0.25±0.03 


0.53±0.07 


0.63±0.08 


0.55 


0.75 


0.10 


0.15 


0.26±0.09 


0.43±0.13 


0.49±0.15 


0.56±0.19 






0.15 


0.20 


0.21±0.06 


0.53±0.08 


0.70±0.08 


0.92±0.09 






0.20 


0.25 


0.43±0.07 


59±0 06 


79±0 06 


1.02±0.09 






0.25 


0.30 


0.21±0.04 


0.62±0.05 


0.84±0.07 


0.98±0.09 






0.30 


0.35 


0.22±0.03 


0.52±0.04 


0.78±0.04 


0.86±0.04 






0.35 


0.40 


0.25±0.04 


0.44±0.03 


0.67±0.03 


0.81±0.05 






0.40 


0.45 


0.18±0.03 


0.43±0.03 


0.65±0.04 


0.81±0.04 






0.45 


0.50 


0.13±0.02 


0.39±0.03 


0.60±0.03 


0.79±0.04 






0.50 


0.60 


0.14±0.02 


0.28±0.02 


0.54±0.03 


0.65±0.05 






0.60 


0.70 


0.09±0.02 


0.22±0.02 


0.44±0.04 


0.52±0.04 






0.70 


0.80 


0.06±0.02 


0.18±0.02 


0.34±0.05 


0.46±0.05 


0.75 


0.95 


0.10 


0.15 


0.24±0.07 


47±0 10 

v / • i j — i — v / ■ iv/ 


54±0 13 


0.55±0.15 






0.15 


0.20 


0.21±0.04 


0.63±0.06 


0.83±0.07 


0.88±0.07 






0.20 


0.25 


0.21±0.04 


0.52±0.04 


0.79±0.05 


0.93±0.07 






0.25 


0.30 


0.21±0.04 


0.49±0.03 


0.70±0.04 


0.85±0.06 






0.30 


0.35 


0.23±0.04 


0.40±0.03 


0.65±0.04 


0.83±0.06 






0.35 


0.40 


0.23±0.03 


0.34±0.02 


0.60±0.03 


0.67±0.04 






0.40 


0.45 


0.23±0.03 


0.30±0.02 


0.48±0.03 


0.54±0.03 






0.45 


0.50 


0.15±0.03 


0.28±0.02 


0.41±0.02 


0.48±0.03 






0.50 


0.60 


0.09±0.02 


0.25±0.02 


0.32±0.02 


0.40±0.03 






0.60 


0.70 


0.06±0.01 


0.17±0.02 


0.26±0.02 


0.31±0.03 


0.95 


1.15 


0.10 


0.15 


0.30±0.07 


0.53±0.10 


0.65±0.11 


0.69±0.14 






0.15 


0.20 


0.35±0.05 


0.65±0.06 


0.80±0.06 


0.86±0.07 






0.20 


0.25 


0.31±0.04 


0.49±0.04 


0.74±0.04 


0.74±0.06 






0.25 


0.30 


0.21±0.03 


0.44±0.03 


0.61±0.03 


0.73±0.05 






0.30 


0.35 


0.13±0.02 


0.36±0.03 


0.50±0.03 


0.57±0.04 






0.35 


0.40 


0.11±0.02 


0.25±0.02 


0.44±0.02 


0.48±0.03 






0.40 


0.45 


0.10±0.02 


0.20±0.01 


0.35±0.02 


0.43±0.03 






0.45 


0.50 


0.10±0.02 


0.17±0.01 


0.27±0.02 


0.37±0.02 






0.50 


0.60 


0.08±0.02 


0.13±0.01 


0.20±0.01 


0.25±0.02 
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^min 




-Pmin 


Pmax 




d 2 a^ 


/(dpd8) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c rad)) 












3 GeV/c 


5 GeV/c 


8 GcV/c 


12 GeV/c 


1.15 


1.35 


0.10 


0.15 


0.29±0.07 


0.60±0.11 


0.68±0.13 


0.65±0.15 






0.15 


0.20 


0.36±0.05 


0.57±0.05 


0.70±0.06 


0.91±0.09 






0.20 


0.25 


0.27±0.04 


0.47±0.04 


0.66±0.04 


0.74±0.05 






0.25 


0.30 


0.25±0.04 


0.42±0.03 


0.51±0.03 


0.59±0.04 






0.30 


0.35 


0.20±0.03 


0.29±0.03 


0.37±0.03 


0.46±0.03 






0.35 


0.40 


0.11±0.02 


0.20±0.02 


0.30±0.02 


0.38±0.03 






0.40 


0.45 


0.05±0.01 


0.16±0.01 


0.24±0.02 


0.28±0.02 






0.45 


0.50 


0.03±0.01 


0.12±0.01 


0.18±0.01 


0.20±0.02 


1.35 


1.55 


0.10 


0.15 


0.33±0.10 


0.56±0.12 


0.66±0.15 


0.73±0.18 






0.15 


0.20 


0.36±0.05 


0.49±0.05 


0.68±0.06 


0.85±0.08 






0.20 


0.25 


0.25±0.04 


0.35±0.03 


0.54±0.04 


0.68±0.05 






0.25 


0.30 


0.21±0.03 


0.28±0.03 


0.40±0.03 


0.47±0.04 






0.30 


0.35 


0.17±0.03 


0.18±0.02 


0.30±0.02 


0.32±0.03 






0.35 


0.40 


0.08±0.02 


0.14±0.01 


0.21±0.02 


0.21±0.02 






0.40 


0.45 


0.04±0.01 


0.09±0.01 


0.16±0.01 


0.16±0.01 






0.45 


0.50 


0.02±0.01 


0.06±0.01 


0.13±0.01 


0.13±0.01 


1.55 


1.75 


0.10 


0.15 


0.25±0.07 


0.57±0.12 


0.66±0.16 


0.78±0.20 






0.15 


0.20 


0.32±0.05 


0.43±0.05 


0.61±0.05 


0.70±0.06 






0.20 


0.25 


0.23±0.04 


0.31±0.03 


0.41±0.03 


0.47±0.04 






0.25 


0.30 


0.11±0.03 


0.20±0.02 


0.30±0.02 


0.35±0.03 






0.30 


0.35 


0.04±0.01 


0.13±0.01 


0.21±0.02 


0.24±0.02 






0.35 


0.40 


0.03±0.01 


0.10±0.01 


0.15±0.01 


0.16±0.01 






0.40 


0.45 


0.03±0.01 


0.07±0.01 


0.11±0.01 


0.10±0.01 






0.45 


0.50 


0.01±0.01 


0.05±0.01 


0.08±0.01 


0.07±0.01 


1.75 


1.95 


0.10 


0.15 


0.19±0.05 


0.52±0.09 


0.55±0.10 


0.68±0.12 






0.15 


0.20 


0.26±0.04 


0.43±0.03 


0.51±0.03 


0.58±0.04 






0.20 


0.25 


0.19±0.03 


0.26±0.02 


0.31±0.02 


0.37±0.03 






0.25 


0.30 


0.08±0.02 


0.16±0.02 


0.20±0.02 


0.24±0.02 






0.30 


0.35 


0.02±0.01 


0.10±0.01 


0.14±0.01 


0.15±0.02 






0.35 


0.40 


0.03±0.01 


0.06±0.01 


0.10±0.01 


0.10±0.01 






0.40 


0.45 


0.01±0.01 


0.04±0.01 


0.07±0.01 


0.07±0.01 






0.45 


0.50 




0.02±0.01 


0.05±0.01 


0.06±0.01 


1.95 


2.15 


0.10 


0.15 


0.25±0.06 


0.40±0.06 


0.47±0.07 


0.52±0.08 






0.15 


0.20 


0.13±0.03 


0.28±0.03 


0.43±0.03 


0.47±0.04 






0.20 


0.25 


0.12±0.03 


0.14±0.01 


0.28±0.02 


0.29±0.03 






0.25 


0.30 


0.04±0.02 


0.10±0.01 


0.16±0.02 


0.12±0.02 






0.30 


0.35 


O.OliO.Ol 


0.07±0.01 


0.08±0.01 


0.11±0.01 






0.35 


0.40 


0.02±0.01 


0.05±0.01 


0.04±0.01 


0.07±0.01 






0.40 


0.45 


0.03±0.02 


0.03±0.01 


0.03±0.01 


0.06±0.01 






0.45 


0.50 




0.02±0.01 


0.02±0.01 


0.05±0.01 
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Table 8: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 (j n /{dpdO) for tin. Each row refers to a different (p mul < p < Pmax, #min — ^ < Vax) bin, where p 
and 9 are the pion momentum and polar angle, respectively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are given. 





"max 




Pmax 




d 2 o* + 


/{dpdff) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GcV/c rad)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.35 


0.55 


0.10 


0.15 


0.07±0.08 


0.53±0.22 


1 06±0 40 

-LtUU — 1 — \J ■ IV' 


1.55±0.49 






0.15 


0.20 


0.29±0.12 


63±0 14 


1.15±0.17 


1.41±0.21 






0.20 


0.25 


0.46±0.12 


0.80±0.09 


1.29±0.12 


1.57±0.17 






0.25 


0.30 


0.65±0.12 


1.02±0.10 


1.44±0.12 


1.90±0.14 






0.30 


0.35 


0.53±0.09 


0.92±0.07 


1.54±0.11 


1.99±0.14 






0.35 


0.40 


0.50±0.09 


1.01±0.09 


1.66±0.11 


1.86±0.10 






0.40 


0.45 


0.43±0.07 


1.07±0.07 


1.53±0.09 


2.07±0.20 






0.45 


0.50 


0.42±0.06 


1.02±0.06 


1.43±0.07 


1.92±0.11 






0.50 


0.60 


0.47±0.06 


0.87±0.06 


1.40±0.08 


1.92±0.12 






0.60 


0.70 


0.36±0.06 


0.74±0.08 


1.22±0.13 


1.77±0.18 






0.70 


0.80 


0.26±0.06 


0.42±0.07 


0.89±0.14 


1.32±0.19 


0.55 


0.75 


0.10 


0.15 


0.23±0.13 


0.65±0.17 


0.87±0.21 


1.19±0.27 






0.15 


0.20 


0.50±0.10 


1.05±0.10 


1.43±0.10 


1.59±0.15 






0.20 


0.25 


0.60±0.10 


1 15±0 09 


1.62±0.15 


2.01±0.18 






0.25 


0.30 


0.68±0.08 


1.05±0.08 


1.58±0.09 


1.99±0.16 






0.30 


0.35 


0.50±0.06 


1.09±0.08 


1.48±0.08 


1.87±0.12 






0.35 


0.40 


0.49±0.07 


0.96±0.07 


1.47±0.08 


1.84±0.11 






0.40 


0.45 


0.46±0.06 


0.86±0.05 


1.29±0.06 


2.01±0.11 






0.45 


0.50 


0.39±0.06 


0.75±0.05 


1.21±0.06 


1.70±0.11 






0.50 


0.60 


0.27±0.04 


0.57±0.05 


0.97±0.06 


1.36±0.10 






0.60 


0.70 


0.18±0.03 


0.37±0.04 


0.70±0.09 


1.00±0.11 






0.70 


0.80 


0.12±0.03 


0.24±0.05 


0.43±0.08 


0.58±0.11 


0.75 


0.95 


0.10 


0.15 


0.50±0.11 


86±0 14 

v / ■ l ' \j — i — v ) ■ i i 


1 03±0 16 


1.14±0.19 






0.15 


0.20 


0.67±0.10 


1.12±0.09 


1.60±0.10 


1.88±0.14 






0.20 


0.25 


0.65±0.08 


1.24±0.09 


1.55±0.10 


1.86±0.13 






0.25 


0.30 


0.60±0.07 


1.02±0.06 


1.34±0.09 


1.77±0.13 






0.30 


0.35 


0.52±0.06 


0.89±0.06 


1.33±0.07 


1.65±0.09 






0.35 


0.40 


0.51±0.06 


0.74±0.05 


1.18±0.07 


1.29±0.06 






0.40 


0.45 


0.34±0.05 


0.57±0.04 


0.95±0.05 


1.16±0.07 






0.45 


0.50 


0.27±0.03 


0.47±0.03 


0.77±0.05 


1.01±0.06 






0.50 


0.60 


0.20±0.03 


0.33±0.03 


0.55±0.05 


0.79±0.06 






0.60 


0.70 


0.09±0.02 


0.19±0.03 


0.31±0.05 


0.49±0.07 


0.95 


1.15 


0.10 


0.15 


0.46±0.09 


0.85±0.11 


1.04±0.14 


1.21±0.17 






0.15 


0.20 


0.63±0.09 


1.01±0.09 


1.45±0.09 


1.83±0.12 






0.20 


0.25 


0.56±0.07 


1.11±0.07 


1.39±0.10 


1.58±0.09 






0.25 


0.30 


0.38±0.05 


0.81±0.06 


1.16±0.06 


1.50±0.11 






0.30 


0.35 


0.24±0.04 


0.55±0.05 


0.90±0.07 


1.18±0.08 






0.35 


0.40 


0.21±0.04 


0.47±0.04 


0.76±0.05 


0.94±0.06 






0.40 


0.45 


0.16±0.03 


0.43±0.03 


0.67±0.04 


0.81±0.05 






0.45 


0.50 


0.12±0.03 


0.32±0.03 


0.48±0.04 


0.64±0.05 






0.50 


0.60 


0.06±0.02 


0.20±0.03 


0.28±0.03 


0.43±0.05 
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^min 
(rad) 


(rad) 


-Pmin 
(GeV/c) 


Pmax 
(GeV/c) 




d 2 a* + /{dpd6) 
(barn/(GeV/c rad)) 












3 GeV/c 


5 GeV/c 


8 GcV/c 


12 GeV/c 


1.15 


1.35 


0.10 


0.15 


0.58±0.12 


0.78±0.11 


1.00±0.13 


1.17±0.18 






0.15 


0.20 


0.74±0.09 


1.05±0.08 


1.34±0.10 


1.53±0.13 






0.20 


0.25 


0.51±0.08 


0.81±0.06 


1.15±0.07 


1.50±0.09 






0.25 


0.30 


0.43±0.07 


0.63±0.05 


0.87±0.05 


1.03±0.07 






0.30 


0.35 


0.26±0.04 


0.47±0.04 


0.63±0.04 


0.81±0.06 






0.35 


0.40 


0.16±0.03 


0.33±0.03 


0.46±0.03 


0.60±0.04 






0.40 


0.45 


0.11±0.02 


0.24±0.02 


0.35±0.02 


0.47±0.03 






0.45 


0.50 


0.09±0.02 


0.19±0.02 


0.26±0.03 


0.35±0.03 


1.35 


1.55 


0.10 


0.15 


0.61±0.14 


0.79±0.14 


1.18±0.20 


1.38±0.24 






0.15 


0.20 


0.76±0.08 


1.12±0.10 


1.41±0.09 


1.66±0.15 






0.20 


0.25 


0.47±0.06 


0.86±0.06 


1.06±0.08 


1.35±0.09 






0.25 


0.30 


0.30±0.05 


0.51±0.05 


0.73±0.05 


0.94±0.08 






0.30 


0.35 


0.15±0.03 


0.37±0.03 


0.55±0.04 


0.58±0.05 






0.35 


0.40 


0.09±0.02 


0.25±0.03 


0.36±0.03 


0.47±0.04 






0.40 


0.45 


0.06±0.01 


0.14±0.02 


0.22±0.02 


0.39±0.04 






0.45 


0.50 


0.03±0.01 


0.09±0.01 


0.15±0.02 


0.23±0.03 


1.55 


1.75 


0.10 


0.15 


0.73±0.13 


0.86±0.15 


1.13±0.19 


1.20±0.22 






0.15 


0.20 


0.68±0.08 


1.02±0.09 


1.34±0.09 


1.52±0.12 






0.20 


0.25 


0.43±0.06 


0.70±0.06 


0.91±0.06 


1.03±0.08 






0.25 


0.30 


0.18±0.04 


0.42±0.05 


0.56±0.05 


0.68±0.06 






0.30 


0.35 


0.09±0.02 


0.24±0.02 


0.39±0.03 


0.43±0.04 






0.35 


0.40 


0.06±0.02 


0.20±0.03 


0.24±0.02 


0.28±0.03 






0.40 


0.45 


0.04±0.02 


0.10±0.02 


0.15±0.02 


0.18±0.02 






0.45 


0.50 


0.04±0.01 


0.06±0.01 


0.09±0.02 


0.10±0.02 


1.75 


1.95 


0.10 


0.15 


0.69±0.11 


0.86±0.11 


0.91±0.12 


0.91±0.12 






0.15 


0.20 


0.49±0.06 


0.74±0.06 


1.02±0.05 


1.02±0.07 






0.20 


0.25 


0.35±0.05 


0.46±0.04 


0.63±0.05 


0.74±0.05 






0.25 


0.30 


0.12±0.04 


0.22±0.03 


0.34±0.03 


0.42±0.05 






0.30 


0.35 


0.03±0.02 


0.13±0.02 


0.22±0.02 


0.22±0.02 






0.35 


0.40 


0.01±0.01 


0.08±0.01 


0.14±0.02 


0.13±0.02 






0.40 


0.45 




0.05±0.01 


0.07±0.01 


0.08±0.01 






0.45 


0.50 




0.02±0.01 


0.04±0.01 


0.05±0.01 


1.95 


2.15 


0.10 


0.15 


0.47±0.09 


0.64±0.08 


0.69±0.07 


0.69±0.08 






0.15 


0.20 


0.44±0.07 


0.64±0.05 


0.66±0.04 


0.80±0.06 






0.20 


0.25 


0.27±0.05 


0.38±0.04 


0.40±0.03 


0.50±0.05 






0.25 


0.30 


0.08±0.03 


0.18±0.03 


0.22±0.02 


0.25±0.03 






0.30 


0.35 


0.02±0.01 


0.10±0.01 


0.12±0.01 


0.15±0.02 






0.35 


0.40 


O.OliO.Ol 


0.05±0.01 


0.09±0.01 


0.08±0.02 






0.40 


0.45 


0.01±0.01 


0.03±0.01 


0.05±0.01 


0.05±0.01 






0.45 


0.50 




0.02±0.01 


0.02±0.01 


0.03±0.01 
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Table 9: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPcr n /(dpdO) for tin. Each row refers to a different (p m i n < p < Pmax, "min — ^ < ^max) bin, where p 
and 9 are the pion momentum and polar angle, respectively. The central value as well as the square-root 
of the diagonal elements of the covariance matrix are given. 





"max 




Pmax 




d 2 a^ 


/{dpdO) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c rad)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.35 


0.55 


0.10 


0.15 


0.28±0.20 


40±0 20 


91±0 36 


1.16±0.48 






0.15 


0.20 


0.30±0.13 


76±0 14 

V / ■ J \J — 1 — \J ■ 1 1 


1.21±0.20 


1.56±0.25 






0.20 


0.25 


0.28±0.09 


0.70±0.11 


1.31±0.12 


1.77±0.20 






0.25 


0.30 


0.46±0.10 


0.76±0.08 


1.38±0.10 


1.85±0.14 






0.30 


0.35 


0.32±0.06 


0.72±0.06 


1.26±0.07 


1.78±0.10 






0.35 


0.40 


0.20±0.04 


0.64±0.06 


1.22±0.07 


1.56±0.08 






0.40 


0.45 


0.29±0.07 


0.62±0.04 


1.08±0.06 


1.37±0.07 






0.45 


0.50 


0.34±0.06 


0.55±0.04 


1.05±0.05 


1.32±0.07 






0.50 


0.60 


0.17±0.04 


0.52±0.04 


0.90±0.05 


1.26±0.08 






0.60 


0.70 


0.17±0.03 


0.43±0.05 


0.84±0.06 


1.26±0.10 






0.70 


0.80 


0.16±0.04 


0.31±0.04 


0.66±0.08 


0.98±0.14 


0.55 


0.75 


0.10 


0.15 


0.28±0.13 


0.50±0.15 


0.96±0.22 


1.15±0.30 






0.15 


0.20 


0.37±0.08 


l.OOiO.ll 


1.41±0.11 


1.58±0.21 






0.20 


0.25 


0.37±0.08 


83±0 08 


1 26±0 08 


1.80±0.14 






0.25 


0.30 


0.43±0.08 


0.80±0.06 


1.27±0.08 


1.52±0.10 






0.30 


0.35 


0.30±0.06 


0.63±0.04 


1.26±0.08 


1.60±0.12 






0.35 


0.40 


0.32±0.06 


0.61±0.04 


1.08±0.05 


1.53±0.09 






0.40 


0.45 


0.23±0.05 


0.60±0.04 


0.97±0.04 


1.27±0.08 






0.45 


0.50 


0.18±0.03 


0.54±0.04 


0.89±0.04 


1.14±0.06 






0.50 


0.60 


0.18±0.03 


0.49±0.03 


0.79±0.04 


1.04±0.06 






0.60 


0.70 


0.15±0.03 


0.39±0.05 


0.59±0.05 


0.80±0.08 






0.70 


0.80 


0.13±0.04 


0.25±0.04 


0.46±0.06 


0.61±0.08 


0.75 


0.95 


0.10 


0.15 


0.31±0.09 


76±0 12 


98±0 14 


1.26±0.19 






0.15 


0.20 


0.49±0.08 


0.91±0.07 


1.42±0.10 


1.76±0.14 






0.20 


0.25 


0.20±0.04 


0.73±0.06 


1.30±0.08 


1.66±0.13 






0.25 


0.30 


0.33±0.06 


0.75±0.06 


1.16±0.06 


1.41±0.09 






0.30 


0.35 


0.32±0.05 


0.69±0.05 


0.96±0.05 


1.24±0.07 






0.35 


0.40 


0.26±0.04 


0.53±0.04 


0.89±0.05 


1.08±0.06 






0.40 


0.45 


0.16±0.03 


0.46±0.03 


0.77±0.04 


0.89±0.05 






0.45 


0.50 


0.13±0.02 


0.38±0.03 


0.60±0.03 


0.88±0.06 






0.50 


0.60 


0.13±0.03 


0.31±0.02 


0.50±0.03 


0.76±0.05 






0.60 


0.70 


0.09±0.02 


0.23±0.03 


0.37±0.04 


0.52±0.07 


0.95 


1.15 


0.10 


0.15 


0.36±0.08 


0.64±0.10 


0.86±0.11 


1.28±0.15 






0.15 


0.20 


0.56±0.10 


0.87±0.06 


1.25±0.08 


1.65±0.12 






0.20 


0.25 


0.46±0.06 


0.78±0.06 


1.07±0.06 


1.39±0.08 






0.25 


0.30 


0.36±0.05 


0.69±0.05 


0.99±0.06 


1.20±0.07 






0.30 


0.35 


0.22±0.04 


0.55±0.04 


0.82±0.04 


0.95±0.06 






0.35 


0.40 


0.14±0.03 


0.38±0.03 


0.68±0.04 


0.77±0.05 






0.40 


0.45 


0.10±0.02 


0.28±0.02 


0.52±0.03 


0.65±0.04 






0.45 


0.50 


0.09±0.02 


0.27±0.02 


0.42±0.03 


0.54±0.03 






0.50 


0.60 


0.08±0.02 


0.21±0.02 


0.30±0.02 


0.41±0.03 
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^min 




-Pmin 


Pmax 




d 2 a^ 


/(dpd8) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c rad)) 












3 GeV/c 


5 GcV/c 


8 GcV/c 


12 GeV/c 


1.15 


1.35 


0.10 


0.15 


0.43±0.09 


0.58±0.09 


0.91±0.13 


1.15±0.19 






0.15 


0.20 


0.46±0.06 


0.90±0.07 


1.21±0.07 


1.49±0.11 






0.20 


0.25 


0.37±0.06 


0.76±0.05 


0.98±0.05 


1.23±0.08 






0.25 


0.30 


0.32±0.06 


0.54±0.04 


0.86±0.05 


1.03±0.07 






0.30 


0.35 


0.17±0.03 


0.36±0.03 


0.64±0.05 


0.80±0.07 






0.35 


0.40 


0.18±0.04 


0.30±0.02 


0.44±0.03 


0.60±0.04 






0.40 


0.45 


0.10±0.03 


0.25±0.02 


0.34±0.02 


0.48±0.03 






0.45 


0.50 


0.06±0.02 


0.18±0.02 


0.26±0.02 


0.36±0.03 


1.35 


1.55 


0.10 


0.15 


0.34±0.08 


0.72±0.12 


1.05±0.17 


1.47±0.27 






0.15 


0.20 


0.52±0.07 


0.99±0.08 


1.27±0.09 


1.63±0.12 






0.20 


0.25 


0.40±0.06 


0.72±0.06 


0.84±0.06 


1.17±0.09 






0.25 


0.30 


0.28±0.05 


0.47±0.05 


0.69±0.05 


0.89±0.08 






0.30 


0.35 


0.13±0.03 


0.28±0.03 


0.52±0.05 


0.59±0.05 






0.35 


0.40 


0.09±0.02 


0.22±0.02 


0.33±0.02 


0.40±0.03 






0.40 


0.45 


0.07±0.02 


0.16±0.02 


0.26±0.02 


0.30±0.03 






0.45 


0.50 


0.05±0.02 


0.12±0.01 


0.19±0.02 


0.22±0.02 


1.55 


1.75 


0.10 


0.15 


0.38±0.08 


0.63±0.11 


0.88±0.14 


1.33±0.24 






0.15 


0.20 


0.43±0.06 


0.85±0.07 


1.23±0.08 


1.37±0.11 






0.20 


0.25 


0.21±0.04 


0.57±0.05 


0.70±0.05 


0.85±0.07 






0.25 


0.30 


0.14±0.03 


0.34±0.03 


0.48±0.04 


0.59±0.05 






0.30 


0.35 


0.10±0.03 


0.25±0.03 


0.36±0.03 


0.47±0.04 






0.35 


0.40 


0.05±0.02 


0.16±0.02 


0.23±0.02 


0.34±0.03 






0.40 


0.45 


0.03±0.01 


0.10±0.01 


0.16±0.01 


0.25±0.03 






0.45 


0.50 


0.02±0.01 


0.08±0.01 


0.11±0.01 


0.18±0.02 


1.75 


1.95 


0.10 


0.15 


0.38±0.07 


0.58±0.07 


0.79±0.09 


0.92±0.13 






0.15 


0.20 


0.38±0.06 


0.60±0.05 


0.83±0.05 


1.03±0.07 






0.20 


0.25 


0.17±0.04 


0.36±0.03 


0.54±0.03 


0.56±0.05 






0.25 


0.30 


0.14±0.04 


0.21±0.02 


0.37±0.03 


0.39±0.04 






0.30 


0.35 


0.11±0.03 


0.14±0.02 


0.25±0.02 


0.25±0.03 






0.35 


0.40 


0.06±0.02 


0.09±0.01 


0.18±0.02 


0.17±0.02 






0.40 


0.45 


0.03±0.02 


0.08±0.01 


0.12±0.01 


0.15±0.02 






0.45 


0.50 


0.02±0.01 


0.06±0.01 


0.08±0.01 


0.12±0.01 


1.95 


2.15 


0.10 


0.15 


0.32±0.06 


0.58±0.07 


0.69±0.06 


0.78±0.08 






0.15 


0.20 


0.22±0.05 


0.53±0.05 


0.66±0.04 


0.77±0.06 






0.20 


0.25 


0.10±0.03 


0.26±0.03 


0.39±0.03 


0.46±0.04 






0.25 


0.30 


0.12±0.04 


0.17±0.02 


0.26±0.03 


0.34±0.03 






0.30 


0.35 


0.05±0.03 


0.10±0.02 


0.16±0.01 


0.19±0.03 






0.35 


0.40 


O.OliO.Ol 


0.06±0.01 


0.11±0.01 


0.10±0.01 






0.40 


0.45 


O.OliO.Ol 


0.04±0.01 


0.08±0.01 


0.08±0.01 






0.45 


0.50 


O.OliO.Ol 


0.04±0.01 


0.05±0.01 


0.05±0.01 
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